Time-evolution of excitations in normal Fermi liquids 
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We inspect the initial and the long time evolution of excitations a Fermi liquids by analyzing the time behavior 
of the electron spectral function. Focusing on the short-time limit we study the electron-boson model for the 
homogenous electron gas and apply the first order (in boson propagator) cumulant expansion of the electron 
Green's function. In addition to a quadratic decay in time upon triggering the excitation, we identify non- 
analytic terms in the time expansion similar to those found in the Fermi edge singularity phenomenon. We also 
demonstrate that the exponential decay in time in the long-time limit is inconsistent with the GW approximation 
for the self-energy. The background for this is the Paley-Wiener theorem of complex analysis. To reconcile 
with the Fermi liquid behavior an inclusion of higher order diagrams (in the screened Coulomb interaction) is 
required. 
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I. INTRODUCTION 

Recently, short-time dynamics of excitations in a Fermi liq- 
uid has gained an increased attention due to the feasibility of 
new spectroscopic techniques [1-4] capable of accessing the 
attosecond time regime. For example in [2] a time delay in 
the range of 80 attoseconds between photoemission from the 
core-level of tungsten and from its conduction band has been 
measured. A further attosecond technique relies on the initial 
excitation of the system by an attosecond pulse and then trac- 
ing the excitation evolution by monitoring the response to a 
second phase-locked laser pulse [4]. This delivers a view on 
how quasiparticle states develop and decay in time on a scale 
well below the time that one may extract from their spectral 
width. The existent experiments on the attosecond time de- 
lay in photoemission call not only for a theoretical determina- 
tion of the measured time delays but most importantly pose a 
question of what new physics can be gained from these highly 
sophisticated experiments. Does the spectral width of the QP 
peak encompass the complete information on how the QP is 
born/decays? As shown below, this is indeed not the case. The 
initial stage of the QP evolution follows a different time law 
as in the long time limit and the decay constants have a spe- 
cific and materials dependent nature, endorsing thus the novel 
physics that can be gained with the attosecond metrology. 

Theoretically, several methods can be used to explore the 
system's dynamics numerically: the density matrix approach, 
the time-dependent density functional or density matrix renor- 
malization group theories and the non-equilibrium Green's 
functions [5-10]. For the latter approach, which can be writ- 
ten in the form the Kadanoff-Baym equations, the two-times 
lesser (greater) Green's functions G s (fi,?2) defined on the 
Keldysh time-loop contour are the central quantities. In the 
equilibrium situation these functions are directly related via 
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the so-called Kubo-Martin-Schwinger boundary conditions to 
the advanced and retarded Green's functions multiplied by 
the corresponding electron/hole distribution functions. In the 
nonequilibrium case they are complicated two-times quanti- 
ties which can be found only approximately by using some 
approximations for the electron self-energy [11, 12]. In order 
to judge the accuracy of such approximations it is desirable to 
have an analytic solution for some limiting cases. 

On the basis of spectral moments calculations [13] for a 
3D electron gas we conjectured an interpolative form of the 
electron spectral function at the energy e: 

A(t; e) = A QP (t; e) exp |-y(e) f + ^ j , (1) 

with Aq P {P, e) the oscillatory quasiparticle part. At longer 
times A(f; e) decays exponentially as expected from very gen- 
eral considerations based on the Landau theory of Fermi liq- 
uids and corrects a spurious divergence of the second spectral 
moment. The latter can be traced back to the fact that the ex- 
ponential decay requires a certain time (r(e)) to set in and that 
initially the decay is quadratic: 

-A(t;e)—* -cr 2 (e)t (2) 

The spectral function, being essentially the decisive equilib- 
rium property, also enters the non-equilibrium two-times dy- 
namics via the generalized Kadanoff-Baym Ansatz [14, 15]: 

-iG f (tuti) = G r (f, - f 2 )G s (f 2 ) - G*(»i)G B (f! - t 2 ), (3) 

and, thus, can be used to calibrate the approximate numeri- 
cal solutions. Although being quite natural, our underlying 
assumptions leading to (1) need to be rigorously verified, a 
task tackled here. In addition, it is desirable to quantify the 
parameters appearing in (1) in terms of measurable physical 
quantities and to offer a scheme for their computations. These 
are the goals of the present work. 
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The cumulant expansion is a well established procedure to 
study the dynamics of many-body systems in the time do- 
main. It amounts to writing the electron Green's function in 
the form: 



g(k,t) = g\k,t)e w \ 



(4) 



The approach gained its wide recognition after Nozieres and 
de Dominicis [16] demonstrated an exact solution of a com- 
plex integral equation for the cumulant function C(k, f) for the 
Fermi edge singularity model. Initial applications to model 
systems include the Mahan's treatment [17] of the Frohlich 
Hamiltonian, the Langreth's study of the singularities in the 
X-ray spectra of metals [18], or the fourth-order cumulant 
expansion for the Holstein model by Gunnarsson et al. [19]. 
More recently the approach was also applied, among others, 
to describe phenomena in realistic systems such as multiple 
plasmon satellites in Na and Al spectral functions [20] or in 
the valence photoemission of semiconductors [21]. Also sim- 
ilarity of this method to the coupled-cluster method broadly 
used in quantum chemistry is well known [22]. 

Typically the method is applied to systems which allow a 
distinct separation of the Hamiltonian into the parts allow- 
ing for the analytical treatment and a coupling that needs to 
be treated perturbatively. A generic example is provided by 
the electron-boson Hamiltonian describing a fermionic sub- 
system (the usual quasiparticles) interacting with the bosonic 
excitations (e.g. phonon or the density fluctuations as will be 
considered below): 



(5) 



k,k' q 



where we narrowed the domain of quantum numbers charac- 
terizing the system to a single wave-vector k as in the case 
of the homogeneous electron gas model. However, the for- 
malism can easily be extended to realistic systems. Here a) q 
describes the energies of bosonic excitations, e& = k 2 /2 is the 
usual particle dispersion in a weakly interacting Fermi liquid 
and V q kk , is the coupling potential. We will use the concept 
of long-lived fermionic excitations (quasiparticles) as a defin- 
ing property of the normal Fermi liquid state [23]. In fact this 
requirement is quite restrictive as it breaks down in, e.g., low- 
dimensional systems. 

Under some circumstances the model (5) is exact: typ- 
ically this is the case when certain matrix elements of the 
Coulomb interactions between a test particle (such as a deep 
core hole [18] or a high-energy photoelectron [24]) are vanish- 
ingly small. For a more general scenario, e.g. as we consider 
here, the accurateness of (5) is less obvious [21]. In view 
of this fact it is interesting to consider the connections with 
other theories. Parallels between the cumulant expansion and 
the many-body perturbation theories (MBPT) in terms of the 
electron self-energy were explored by Aryasetiawan [20] in 
the lowest order. But does this correspondence hold at an ar- 
bitrary order? Development of the formalism and an answer 
to this question will be provided in Sec. II. 

The cumulant function C(k, f), we show, can be written very 
accurately in terms of the dynamical structure factor S(k, to). 



This allows us to obtain exact analytical results for the pref- 
actor cr 2 (e) of the quadratic decay (2). It is interesting that 
also terms proportional to {—iff can be written in a concise 
analytical form. These results do not require any further ap- 
proximations in addition to the ones discussed in Sec. II and 
follow from exactly known sum-rules and the asymptotic be- 
havior of S(k, oj) (Sec. III). 

The quasiparticle uncertainty cr 2 (e) (2) can alternatively be 
obtained from the zeroth spectral moment of the electron self- 
energy [25]. While this approach is perfectly justified for fi- 
nite systems [26] a careful analysis must be done in the case 
of the homogenous electron gas (HEG) model. Here the diffi- 
culties arise from a particular asymptotic behavior of the elec- 
tron self-energy: decay as nT 3 ' 2 for to — > oo and its vanishing 
imaginary part for oj < co*(k). In Sec. IV we trace the ori- 
gins of these features. On the basis of the Paley- Wiener (PW) 
theorem [27] we demonstrate that such restriction of the spec- 
trum implies an unphysical spectral function which in the time 
domain asymptotically decays faster than the exponent. This 
is, however, a consequence of GW approximation for H(k, at). 
The paradox is further resolved by considering higher order 
contributions to the electron self-energy. 



II. CUMULANT VS. SELF-ENERGY EXPANSION 

For the Hamiltonian (5) the interaction between fermionic 
and bosonic degrees of freedom is given by the fluctuation 
potential [24] with matrix elements: 



v k,k' 



An 



k-k' 



|2 ^k-k' 



(6) 



where q is a quantum number characterizing the excited 
bosonic states and p\ is the kth Fourier component of the 
fluctuation density operator between the ground state and a 
state with one boson with the quantum number q excited. The 
choice of the interaction form is not arbitrary: it guarantees 
that the lowest order diagram for @(k, t) in the model (5) cor- 
responds to the GW approximation (for more details on this 
approximation see Sec. IV) for the initial fermionic system. 

It is instructive to derive C(k, f) starting from MBPT and 
using the method of Aryasetiawan [20]. By comparing the ex- 
pansion of the exponential in (4) to that obtained by iterating 
the Dyson's equation: 

0(k,t)=g p (k,t) 

+ jdTjdT'g°(k,t-T)I.(k,T-T')g(k,T') (7) 

we can easily verify that Cik, t) and l,(k, f) should have the 
same lowest order expression in terms of the interaction. 
We separately consider the particle (k > k F ) and the hole 
(k < k F ) cases. The non-interacting Green's function is given 
by g° + (k, i) = -i8(t) e~ iekt and Q a _ (k, t) = iO(-t) e~ iet ' ', respec- 
tively. We represent the full Green's functions as @ ± (k, t) = 
Q° ± (k, t)e c±< -'K These notations are different from the ones 
used by Langreth where they denoted two differently defined 
Green's functions. For the electron self-energy we adopt the 
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standard (non self-consistent) expression (cf. Eq. (25.1) of 
[28]): 

= i f ^(q - M + 6) g°(q, t), 5 -» +0, 
J (27rP 

with the screened Coulomb interaction given by: 

2oj I V 9 ' 2 

'U'U.c.,) = v(k) + y 2 n(fe,w) = r(A) + ^ ' ''' 



w 2 - w 2 



(8) 



where i>(A:) = An/k 2 is the Coulomb potential and Y\(k,oj) is 
the full bosonic propagator or the density -density response 
function in this particular case. The latter is related by the 
fluctuation-dissipation theorem (cf. Eq. 3.74 of [23]) to the 
dynamical structure factor: 

S(k,o>) = — lm Tl(k, oj). 
n 

For the w-dependent part of the screened Coulomb interaction 
we use the spectral representation: 



<W(k,t) = iv 2 (k) d(jjS(k,<jS)e- la ", 
Jo 

where the dynamic structure factor is expressed in terms of 
the imaginary part of the dielectric function (e = s' + is"): 

k 2 s"(k,oj) 
S(k, oj) = — — -J 9(o>). 
An 1 \s(k, oj)\ z 

Clearly, four cases arise depending on the length of q and k. 

a. k > kf and q > kf: Inserting the expressions for £ 
and <W in Eq. (7) we obtain: 



C + (k, t) = - f ^v 2 (q) f dcuS(q - k, oj) e-'" 

Jq>k F (ZX) Jo 

x jj diTT^e^-^-"^-^ 



(9) 



where the integration domain is determined by the condition 
9(t) 0{t-T) 6(t') O(t-t'). This is, in fact, a finite domain which 
can be integrated as follows: 



f(v) = ff d(TT')e iv(T - T,) = f dT f dT'e iv(T 
JJd + Jo Jo 

/ 

Jo 



= dr(t-T)e lVT = 
la 



-T') 

1 + ivt - e m 



(10) 



b. k > k F and q < k F : We have to use the hole prop- 
agator for the intermediate line. This changes the sign of the 
expression and modifies the integration domain which we rep- 
resent as two terms: 

8(t - t) 6(t') 0(t' - r) = 6(t - t) 9(t') (1 - 6»(t - t')) (11) 

The first term here represents an additional contribution perti- 
nent to GW approximation only: 



c. k < kp and q < kf.- The procedure is similar with 
the only difference in the domain of the integration D = 
9(-t) 9(t - t) O(-t') 9(t' - t) which also integrates in terms 
of/(v): 



d(Tr) e 



(V(t-t') 



dr' e 



i'v(t-t') 



: /(V). 



Hence, C_(f) can be written in the same form as (9) and we 
will use C(k, t) as a common symbol for both C ± (k, t). There 
also is a contribution from the intermediate hole line q > k F 
to C_(f) which can be evaluated along the same lines as (12). 
Finally, we redefine the variable for momentum integration as 
q - k — > q and obtain in line with Langreth [18]: 

Jr»oo 
dajS(q,u>) 


x/(e|k| -%+q| -oj,t). (13) 

The central quantity of this study - the dynamical structure 
factor - although expressed almost identically (except for the 
(12) terms) in the many-body perturbation and in the cumulant 
expansion theories, originates from different approximations. 
In the former case it is the vertex function in the expression 
for the self-energy that is neglected, while for the latter it is 
assumed that the Hamiltonian can be written in the electron- 
boson form (5). For the homogenous electron gas model the 
justification mostly comes from MBPT although S(k, a>) can 
be rather complicated function even for simple systems [29]. 



III. SHORT AND LONG TIME LIMITS 

Equation (13) is general enough to treat all the cases pre- 
sented in Tab. I. Compared to Eq. (44) of [18] we additionally 
allow the test particle to scatter (i.e. to exchange its momen- 
tum with the bosonic excitations) whereas a deep core in [18] 
is assumed to have an infinite effective mass. The short-time 
limit of the electron Green's function crucially depends on the 
exact form of fermionic dispersion, on the boundness of the 
bosonic spectrum and on the actual form of the coupling po- 
tential. 

As a first application we compute the leading expansion co- 
efficients of the cumulant function in the short-time limit: 



C(k, t) 



o- 2 (k)_ t2 + £3(fe) f3 + 



2! 



3! 



(14) 



We note, however, that such expansion does not imply analyt- 
icity of the function in vicinity of t = 0. Just the opposite, 
higher expansion coefficients diverge starting from c^ in 2d 
case and from c 7 in 3d based on very general properties of 



dr I dT e lv{T ^ T ) j. (12) the density-density response function (Im H(k, oj) ~ oT^ d ^ 2 , 

00 Jo n 1 39 of T231Y Although where the Hi versrence occurs exactlv 



For the second term in (1 1) we have the same form and sign 
as in (9). Thus, the integration over q < kf can be combined 
with (9) resulting in the sum over all momenta. 



p. 139 of [23]). Although where the divergence occurs exactly 
can be modified by including higher order terms (in S(k, oj)) in 
the expression for the cumulant function, this will not restore 
the analyticity. The prefactor of the quadratic decay (Eq. (2)) 



4 



TABLE I. Electron-boson models and main results for the long and short-time limits of the electron Green's function 



FcIIIllUIllC UlSpeiSlUIl 


uispersiomess pnonons 


L^ispersioniess piasinons 


Th 1 ntrnn rn~*1 *a noire 

j^ieciron-noic pairs 


Deep hole: H = Ec^c 




C(t) ~ e-" "' 
Langreth (1970) 


C(t) ~ -a[ In \Dt\ + in/2 sgnf] 
Langreth (1970) 


Valence states: Hq = 6kc\ck 


C(t) ~ — (2it) 112 e-^ ' 
Mahan (1961) 


C(t) ~ e^"' 
Aryasetiawan (1996) 


this work 



can be computed by evaluating the second derivative of (13) 
at t = 0: 

cr 2 =nJ]v 2 (q)S(q), (15) 

9 

where the static structure factor is defined as 

1 f°° g^Q q 1 

S(q) = - dcoS(g,w) (16) 
n Jo 2<^ P (q) 

It follows then that a 2 is independent of k and coincides with 
the local contribution to the zeroth spectral moment of the 
electron self-energy obtained by Vogt et al. [13]. The long 
wave-length expression in Eq. (16) follows from the exactness 
of the random phase approximation (RPA) in this limit. In the 
opposite case (i.e. q — » oo) the structure factor approaches 
unity, however, the subleading term RPA fails to reproduce. 
In order to accurately compute <x 2 the parameterized structure 
factor of Gori-Giorgi et al. [30] based on the quantum Monte 
Carlo results was used [13]. The convergence of the integral 
(15) is ensured by the limit: 

lim q z+1 [S(q) - 1] = -n2 z ng(0), (17) 

q^oo 

where g(0) is the value of the pair correlation function for two 
electron at the same position and z is the dimensionality of a 
system. 

It is not obvious from the outset that the ailc) coefficient 
should take a finite value: this heavily relies on the exact form 
of the structure factor in the asymptotic (g — » oo) limit. By 
using the /-sum rule: 

dcja>S(q,ti)) = ne q , (18) 

o 

where n is the electron density we obtain: 

c 3 (lc) = -in ^ u 2 (<?)[(e|k| - %+q|) S(q) + e q \. (19) 

q 

In the simplest case of a hole state at the band's bottom the 
convergence of the integral regardless of the system's dimen- 
sion (z) is guaranteed by the limit (17). For k > the term 
linear in k vanishes after the angular integration and we finally 
obtain the ^-independent result: 

c 3 = inJ^u 2 (q)e q [S(q)- l]. 

q 



Finally we notice that the leading terms of Eq. (13) in 
the long time-limit are the constant and the linear ones, i.e., 

C{k, t) - — > y - i €k)t as expected from the exponential 
quasiparticle decay (cf. Eq. (7) of [20]). 

d. Non-analyticity of the spectral function at t — The 
asymptotic behavior of the density-density response function 
at large a> leads to diverging expansion coefficients in Eq. (14). 
This, in turn, gives us a hint that the cumulant function is prob- 
ably nonholomorphic at t — 0. Such property is, however, not 
an exception, but rather the rule as Tab. I demonstrates. We 
will sketch below how all the results presented in this table can 
be obtained in a unified way from Eq. (13) and will also show 
that the same applies to the normal Fermi liquids, in particular 
due to the scattering of valence electrons with the generation 
of electron-hole pairs (viz. "this work" in Tab. I). 

One of the most interesting scenarios is the case of a core 
hole coupled to electron-hole excitations. In the limit of infi- 
nite mass of the fermion and t — > oo there is a singular term 
that arises from the frequency integration in Eq. (13): 

C(f) = -X^7^ du 2 ln|Df| + -sgnf 

Y Hq, 0)| 2 Jo q 2 

where < 77 < \ is called the Anderson singularity index [31] 
which can be given in terms of the scattering phase shifts 61 of 
the statically screened potential CW(q, 0)) as 

In the frequency domain the resulting spectral function ex- 
hibits a singularity ^ which for the finite hole's mass 
and z > 1 is completely washed out by the effect of scatterer 
recoil as was demonstrated by Nozieres [32]. This equiva- 
lently can be seen from our model (13) where the momentum 
angular integration of the function /(%! - f|k+q| -<w , removes 
the singularity. 

The cumulant function resulting from the interaction with 
plasmons has a simple structure which likewise follows from 
(13) by using the limiting form of the structure factor: 

q^Q q 2 <jJ„(0) 2 

S{q,oS) > _ 

8n to p (q) 

In the frequency domain this leads to the main quasiparticle 
peak accompanied by a sequence of satellites displaced by 
no) p (0) [12, 33]. In many realistic materials these are indeed 
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observed features [20, 21]. Recalling our remark in the in- 
troduction concerning the current status of the experimental 
attosecond spectroscopy, it is obvious that these structures are 
interesting candidates for the tracking of the development of 
main and satellite peaks to their static limit. 

Another interesting case, likewise in the t — > oo limit, arises 
from the polar coupling ~ l/q between an electron and the 
non-dispersive optical phonon with the energy a>o [17]. It re- 
sults in the effective structure factor S(q, co) ~ 1 /q 2 6(a> - too). 
There, the non-analytic terms stem from the momentum inte- 
gration: 



1 



(2tt) 3 / 2 



if- 



To see the connection we express asymptotically the spectral 
function as A(k, co) ~ C/w 2+3 ^ 2 and perform the Fourier trans- 
form. Since at co — > -oo the spectral function decays faster, 
in fact on the GW level it is even zero below a certain thresh- 
old value of co, it is sufficient to perform the transform on a 
semi-bounded interval: 



A(k. 



dco C 
2nc~J^ 



Among several resulting terms one has to pick up the one 
independent of the cut-off to c . It exhibits the same time- 
dependence and the density scaling ~ (ar s ) 2 as Eq. (20). 



For the momentum state k = we can use the representation 
of /(e|k| - £jk+q| ~ w , m terms of a single time-integral (10) 
and obtain: 



C(0,0 ~ dr(t-T) . 



e -«"oT „ —(2it) 1/2 e- it00 ' . 

0)Q 



In the opposite case, i. e. t — > 0, the non-analytic terms 
originate from the coupling to the particle-hole (p-h) contin- 
uum. We can split the momentum integration into a finite 
interval q < q c yielding just the well-behaved analytic part 
of C{k, f) and the interval extending to infinity. The value of 
q c can always be chosen large enough so that the real part 
of the dielectric function on the second interval approaches 
unity. This considerably simplifies the dynamical structure 
factor which results now from the imaginary part of the Lind- 
hard formula (Eq. (5.35) of [28] for z = 3) only: 



S(q,co) 



1 1 ar s 
n v(q) q 3 



where the tilde denotes the use of rescaled quantities, i.e. q = 
qk F , co = oj€f, f = t/sf and so on. After the substitution 
co = q 2 + 2qA we can first integrate over the interval |A| < 1. 
This yields a trigonometric expression which is just a constant 
in the lowest order of 2qf: 



I 



doS{q,cj)e- KS >-? )f = -^kl+0(2qf). 



The leading non-analytic term of the remaining momentum 
integral reads: 



(2. 



q 2 dq V \q )e - 2 ^ --U^ 1 , 



Finally we perform the time integration as in (10): 



C(0. 



2ki r> 



T-\l/2 
2lT\ 



rfr(r-r)| T =-^Y[- 



lit 



5/2 



(20) 



Such time-dependence is easy to reconcile with well known 
asymptotic behavior of the electron self-energy as a function 
of frequency [13, 34]: 



~ OJ— >00 

Iml(fc,w) 



16 V2 (ar s ) 2 
3n w 3 / 2 ■ 



IV. SELF-ENERGY IN FREQUENCY SPACE 

In 1965 Lars Hedin formulated a system of functional equa- 
tions [35], carrying by now his name, that relate the elec- 
tron self-energy 11(12), the irreducible polarization propaga- 
tor !P(12), the screened Coulomb interaction ^(12), the ver- 
tex function T(12; 3) and the electron Green's function Q{\2). 
The Hedin equations are becoming one of the major the- 
oretical tools for the treatment of correlated many-particle 
systems [36-38]. The homogeneous electron gas (HEG) in 
two or three dimensions is a prototypic model which allows 
for testing various approximations to the exact Hedin' s equa- 
tions. Earlier applications revealed important features of the 
single -particle spectrum [39^41]. These single-shot calcula- 
tions were extended by several authors to the self-consistent 
level [42, 43], higher order diagrams were included [44, 45]. 

In view of large efforts devoted to the study of these model 
systems it is surprising that some aspects remained unnoticed. 
Thus, it is commonly believed that HEG in two or three di- 
mensions serves as a perfect illustration of the Fermi liquid 
concept [23], that is a many-body fermionic systems with 
long-lived excitations: quasiparticles . Two marked proper- 
ties distinguish them from other excited states: i) they can be 
brought in a direct correspondence with real particles (elec- 
trons) of a fictitious non-interacting many-body system; ii) 
they are characterized by the life-time, which tends to infin- 
ity as the particle's energy approaches the Fermi level (ep). 
It also implies that asymptotically the decay is exponential 
exp(-yf)> with the decay constant being quadratically depen- 
dent on the energy (y(e) ~ e 2 /e^). At e — » £f this constant 
can be computed perturbatively, and it is sufficient to consider 
the lowest-order term giving a non-vanishing imaginary part 
of the self-energy. In view of this it is intriguing that a rigor- 
ous proof can be given that the lowest-order diagram yield the 
spectral function inconsistent with the asymptotic exponential 
decay. 

One recognizes that pronounced features in the spectral 
function appear at energies En that are approximately given 
by Ef, — + S(fe, Ek), where ImZ(A:, co) ~ S(co - e\k+q\ ± ^q)- 
These resonances are surrounded by the incoherent back- 
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FIG. 1. a) First and b) second order hole scattering mechanisms. 
While the two excitations in b) carry in total the same momentum as 
a single excitation in a) the energy transfer is much larger. 



ground which has the same extend as the self-energy: 



A(k,co) = -\lm g(k,to)\ 
n 

1 



\\mY.(k,co)\ 



n | co - e k - Re S(Jk, co)\ 2 + |Im Z(Jt, w)| 2 



(21) 



In the lowest order of the screened interaction a particle can 
only loose its energy (e k ) by generating a sz'ng/e bosonic ex- 
citation co q . Since only a finite momentum can be transfered 
also oj q is finite and, thus, the self-energy has a semi-bounded 
support (limited from below) (Fig. 1). From this, in view of 
(21) follows A(k, co) = for co < of{k). 

The last property allows us to apply the Paley-Wiener (PW) 
theorem [27] which we present below for completeness in the 
original formulation: 

Theorem 1 Let <f>(x) be a real non-negative function not 
equivalent to zero, defined for -oo < x < oo, and of integrable 
square in this range. A necessary and sufficient condition that 
there should exist a real- or complex-valued function F(x) de- 
fined in the same range, vanishing for x > xofor some number 
xq, and such that the Fourier transform G(x) of F(x) should 
satisfy \G(x)\ — <p{x), is that 



f 



i log m\ 

l+x 2 



dx < oo. 



(22) 



The modern formulation [46] only slightly relaxes the condi- 
tions on the functions: 

Theorem 2 For tb{x) e L 2 (R) and tb{x) > the integral (22) 
converges <==> there is a function F £ L 2 (R) with a semi- 
bounded support such that (p = \T[F]\ a. e. in R, and T is the 
Fourier-Plancherel operator. 

From here follows: 



Corollary 1 

F(y) = for y < y 

X — i 

Mx)\ - 



I log 0WI 
1 +X 2 

e\p(-Bx a ), 



dx < oo => 
where a < 1. 



If we identify now F with A(k, co) the deviation from the expo- 
nential decay for A(k, t) follows: a result in a clear contradic- 
tion with the Fermi liquid theory. Before proceeding with the 



resolution of the paradox we present our method for numeri- 
cal calculation of co*(k). It is sufficiently general in the sense 
that there is no limitation on the dimensionality (z) of the sys- 
tem and it is not limited to the first order expression. From 
now on we will only be using rescaled quantitates, tilde will 
be omitted for clarity. In order to make the model amenable 
for the numerics we introduce the following representation of 
the screened Coulomb interaction: 



2 J [co-n(k,A) + 



w(k,A) 



IT] 



oj + Cl(k, A) - ir\ 
(23) 

where w(k,A), Q.(k,A) are some real functions that will be 
specified below. Our representation takes advantage of the 
fact that the imaginary part of the dielectric function and 
the screened Coulomb interaction is different from zero only 
in the stripe area in the cj-k plane and along the plasmonic 
line [23]. The limits for the particle-hole continuum are given 
(for z> 2) by: 

max{0, a>-(k)} < M < 0)+(k), with cj ± (k) = k 2 + 2k, (24) 

where it is convenient to parameterize the trajectories on the 
stripe (24) as: 

a(k,A) = k 2 +2Ak. 

Thus Eq. (23) is nothing but the spectral representation (see 
e. g. Eq. 4 of [42]): 



<W(k, oj) = v(k) + 



f 

Jo 



2o/B(k, to') 



do/, 



The integral over A is to be understood in a generalized sense: 
this parameter can assume both discrete values when we de- 
scribe a single excitation such as plasmon or be a continu- 
ous variable for particle-hole excitations. In the former case 
it reduces to the plasmon model approximation (cf. Eq. 25.1 1 
of [28]): 



nV°(fc, co) = v(k) 



1 + 



u 2 - col{k) 



The bare Coulomb part can also be obtained from (23): con- 
sider the limit Q.(k,A) — > w(k,A) — > oo. The fact that we 
can represent all contributions to W(k, co) in a unified way is 
crucial for our discussion: one does not need to separately 
consider diagrams with bare or renormalized interaction lines. 
The former can be obtained from the general case by formally 
taking the limit of the final expression. 
We write the Green's function as: 



&\k,co) 



+ 



1 - n k 



co - e k -iT] to - e k + ii] 



where n k denotes the occupation of the state with the momen- 
tum k and consider the two lowest order diagrams for the elec- 
tron self-energy OV ]: 

2 (1) (1,2) = ;£ (l,2)OV°(l + ,2), (25a) 
I (2) (l,2) = i 2 jj OV (r,4)£°(l,3)£°(3,4) 

x^°(4,2)nV°(3 + ,2)c/(34). (25b) 
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Our representation of the screened Coulomb interaction al- 
lows to compute the electronic self-energy relatively easy us- 
ing the maple computer algebra system [47]. The final results 
can be recasted in the form of momentum and A integrals over 
complicated domains that we denote as D±, where + desig- 
nates particle (hole) state, and (z) is the order of a diagram. 
Since we are only interested in the phase-space where each 
diagram contributes we skip here the explicit expressions for 
Y(k, oj) and present only the results for 2)® (k, cS): 

D m = n p 5{u) - e p + Q.(q, A)), (26a) 
£) {2a) = (1 - npo ) Hp] n pi 5{u> + e Po - e pi - e Pl ), (26b) 
D {2h) = n po 6{cj - e P0 + «(?!, + «(«2. h)), (26c) 
fl W =n PJ 6(u-e Pl +n(q u A 1 )) 

x[l + fa(quqz,puP2,Po)(l -n Pl ) 

Mqi ,q2,Pi,Pi, Po)(l - n Po )\ + (1 <-> 2). (26d) 

where we introduced the following vectors p = k + q, po = 
k-qi -q2, pi = k-qi, and p2 = k-q2. D_ (k, to) describes 
the simplest first order process when a hole scatters to another 
hole-state hereby generating a plasmon or a p-h pair. D^ 2a) 
describes a hole scattered to a two-holes-one-particle (2h-p) 
state, whereas £) (2fe) stands for a process when a hole looses its 
energy by the generation of two bosonic excitations. The last 
term is not interesting because it just renormalizes D_ (k, to). 
Expressions for the particle states analogous to (26) can be 
obtained by the use of the particle-hole symmetry. 

We compute the integrals involving D ( ±(k, to) by using the 
Monte-Carlo approach and formulas for the momentum in- 
tegration presented in Appendix A (Fig. 2). Thus, the first 
order contribution is obtained by throwing a quartet of ran- 
dom numbers (k, A, q, y) consistent with the integration do- 
main. Henceforth, we verify the condition imposed by the 
^-function and determine possible values of to. As for £p; a) 
the probability distribution is obtained from a set of 5 num- 
bers (k, q, Q), whereas we need to additionally sample 
over A\j, random variables for £) <2i \ 

In agreement with our simple argument we see that the 
phase-space for the first order processes is limited. The same 
is observed in the simplest second order process (it includes 
also contribution from two bare interaction lines) in view of 
the same arguments. The existence of a critical upper momen- 
tum for the plasmons also restricts the phase-space available 
for the h —> h + 2pl scattering. The situation is completely 
different for the /!—>/; + (p-h) + (p-h) events: even though 
the hole can only loose a finite momentum the shares between 
the excitations can be large (Fig lb), resulting in an arbitrarily 
large energy transfer (cf. Eq. (26c)). Hence, the self-energy 
has an unbounded support, the Paley-Wiener theorem cannot 
be applied and the Fermi liquid behavior is restored in the sec- 
ond order. 

Our analysis is also important for practical calculations 
since it allows to determine a priori where a certain diagram 
might contribute. It is interesting to notice a sequence of plas- 
monic peaks in the see of p-h excitations. By expanding the 
cumulant function C(k,t) ~ e la>p ' (Tab. I, third column) and 




k=1.6k F 




FIG. 2. Monte-Carlo calculation for 3d HEG at r s = 5 density of 
the first and second order diagrams contributing to the electron self- 
energy (Eq. (26)). 0^ b) describes 3 second order processes: gener- 
ation of two plasmons, of two particle-hole pairs, or of one plasmon 
and one particle-hole pair. The latter as well as D^ 2c) are not shown 
because they only represent corrections to the first-order processes. 
The Fermi energy is set to zero. 



computing the Fourier transform one sees that their weight 
decays as e~"a" jn\. Guzzo et al. [21] estimated a ~ 0.3 for 
silicon. Therefore, plasmons will only be important at low or- 
ders whereas the tails of the spectral functions are shaped by 
the p-h scattering mechanisms which lead to the power-law 
decay. Where such a crossover occurs depends, of course, on 
the specific system parameters. 

The phase-space arguments provide a partial account of 
the problem. The inclusion of matrix elements can modify 
the self-energy substantially as the comparison of Fig. 2 and 
Fig. 3 shows. This can be best seen at the Fermi level (set to 
zero in our calculations). While both methods lead to a van- 
ishing self-energy in this limit the way how it approaches zero 
is rather different. But how feasible is the realistic calculation 
of next order diagrams? To answer this question let us con- 
sider Eqs. (26). There, the second order terms were evaluated 
by using at most a 7-dimensional sampling. The full-fledged 
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Re Z(co^^ 
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-3-2-10123 

Energy to/to p , 

FIG. 3. Exact first order self-energy of the 3d HEG at r s = 5 den- 
sity. The Fermi energy is set to zero. The real part is shown without 
including the static (exchange) part. The values are e x = -1.490o» p 
and e x = -0.225o> ; , for k = 0.4k F and k = l.6k F , respectively. The 
plasmon energy in the long wave-length limit is a> p = 2.103e f . 



suit from rather inefficient coupling to p-h excitations. We 
also provide a concise analytic form for the expansion coeffi- 
cients of the cumulant function at t — > 0. 

The short-time limit was analyzed using the cumulant ex- 
pansion method, while the asymptotic behavior at longer 
times was studied using the ordinary many-body perturba- 
tion theory. Thus, it was necessary to establish a connec- 
tion between both methods. We have shown that in the first 
order in the screened Coulomb interaction the cumulant ex- 
pansion can be recovered from the MBPT expression for the 
self-energy, although, some terms are additionally present in 
MBPT. For a finite systems we proposed a similar approach 
and presented results supported by full numerical calculations 
for Na-clusters [25] and C6o [26], which supports the general 
nature of the time evolution law of the spectral function. The 
specific, material-dependent, and quantum size effects are en- 
capsulated in the decay constants. 
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Appendix A: Some momentum integrals 



evaluation, in contrast, would require an 8-dimensional inte- 
gration for each k and a> values. In some specific cases simpli- 
fications might be achieved such as in exact analytic treatment 
of the second-order exchange term by Onsager et al. [48]. On 
the other hand, for practical applications some synthetic ap- 
proaches might be promising [29, 49]. 



V. CONCLUSIONS 

In this contribution we performed a detailed analysis of the 
formula A(t; e) = Agp{t; e) exp f— y(e) t+ ^, , j for the time evo- 
lution of the spectral function for extended systems. The vi- 
olations such as i) non-analyticity at short times, ii) reduction 
of the spectral weight of the quasiparticle peak or iii) a non- 
exponential decay at the long-time limit were found. Surpris- 
ingly, our theory reveals that these features are either artifacts 
of approximations used (iii) or are rather weak (i) as they re- 



A single momentum integral leading involving D m {k, oS) 
can be computed as follows: 



dqf(q, |k + q| 2 ) 



n r°° r(k-q) 2 
- qdq 

K J0 J(k+q) 2 



dyf(q,y), (Al) 



with y = |k + q| 2 . The two momenta integrals involing Q y L 
can be computed by introducing symmetrized variables as 
in [48]: 

j d qi j dq 2 /(|k - Ql - q 2 | 2 , |k - qi| 2 , |k - q 2 | 2 ) 

= I I dq I " 2q ' 2 ' lQ " Q ' 2 ' lQ + Q ' 2) 
= — I dq I dyA dQ\ dy 2 Q 

°k Jo J(k+2q) 2 JO JqQ 

x fim ,q 2 + Q 2 ~ 2y 2 , q 2 + Q 2 + 2y 2 ), (A2) 
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where where 

q = k- i(qi +q 2 ), Q=^(qi-q 2 ); 

yi = |k - 2q| 2 , 5/ 2 = (q • Q). 
Analogically for the Tpl b) term we have: 



1 1 

r r q=^(qi+q2), Q=~(qi-q2); 

J d qi J rfq 2 /(|k - qi - q 2 | 2 , |q t | 2 , |q 2 | 2 ) 1 - 



l -^d^dQ /(|k - 2q| 2 , |q + Q| 2 , |q - Q| 2 ) 

n 2 r°° r(k-2q) 2 s*oo , — qQ 

- — I dq I dyi I rfg I ^2 2 

OK Jo J(k+2q) 2 Jo J 9 G 

x , q 2 + Q 2 + 2y 2 , q 2 + Q 2 - 2y 2 ), (A3) 



y x = |k - 2q| 2 , j/ 2 = (q • Q). 
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